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INTRODUCTION 


Huang  et  al.(ref  1)  addressed  the  question,  "Can  the  fractal  dimension  of  images  be 
measured?"  They  found  that  measured  fractal  dimension  D  varied  monotonically  with  true  D  in 
their  mathematically  constructed  fractal  surfaces.  However,  they  found  large  discrepancies 
between  their  measured  and  true  fractal  dimensions  and  concluded  that  the  fractal  dimension  of 
images  could  not  be  measured. 

Recently,  the  effects  of  normally  and  uniformly  distributed  random  noise  on  the  apparent 
fractal  scaling  of  surfaces  in  3-space  was  investigated  by  Meisel  (ref  2).  It  was  found  that  the 
fractal  dimension  of  images  of  the  sort  investigated  in  Reference  1  could  be  measured  and  that 
fractal  analysis  of  surfaces  in  the  presence  of  uniformly  or  normally  distributed  random  noise 
yielded  measured  D  larger  than  the  intrinsic  D  by  predictable  amounts.  The  sensitivity  to  noise 
was  greatest  for  D  near  2.0  and  decreased  smoothly  as  D  approached  3.0.  Quantitative  estimates 
of  the  increase  in  the  measured  D  to  be  expected  for  a  given  amount  of  either  type  of  noise  on 
Brown  surfaces  (refs  3,4)  were  presented.  Thus,  one  could  conclude  that  the  fractal  dimension  of 
images  of  surfaces  in  3-space  in  the  presence  of  normally  and  uniformly  distributed  random  noise 
can  be  measured,  and  estimates  of  the  "errors"  in  D  can  be  made  if  the  magnitude  of  the  noise  is 
known. 


The  problem  of  measurement  of  D  for  images  of  surfaces  in  3-space,  which  are 
contaminated  by  systematic  but  localized  noise,  is  addressed  here.  In  particular,  the  present 
results  suggest  that  one  can  obtain  useful  measures  of  D  by  analysis  of  laser  scanning  confocal 
microscopy  (LSCM)  (ref  5)  images  of  surfaces,  even  when  the  images  are  contaminated  by  a  high 
density  of  "phantom  surface  artifacts." 

Laser  scanning  confocal  microscopy  (ref  5)  is  predicated  on  the  idea  that  in  a  system  with 
a  very  sharp  focus,  the  detected  intensity  will  be  maximal  when  a  piece  of  a  surface  being  imaged 
is  precisely  in  focus.  Thus,  for  example,  in  LSCM  at  fixed  z  =  one  scans  along  the  y  direction 
at  fixed  X  =  Xj  and  records  values  in  an  image  matrix  A  at  y-values,  yj  corresponding  to  "large 
enough"  local  intensity  maxima,  A(xi,  yj)  =  zj^.  After  the  scans  have  been  completed  for  the  entire 
set  of  Xj,  the  value  of  z  is  increased  to  the  next  larger  value  in  the  set  and  the  procedure  is 
repeated. 

A  complication  associated  with  practical  LSCM  surface  characterization  is  the  tendency 
for  images  to  be  contaminated  with  phantom  surface  artifacts.  Phantom  surface  artifacts  appear  in 
spatial  regions  where  no  part  of  a  surface  is  actually  in  focus  in  a  line  scan.  In  such  regions,  a 
practical  LSCM  may  mistake  a  local  intensity  maximum,  which  tends  to  occur  above  or  below  a 
"geometrical  feature"  of  the  true  surface,  as  if  it  were  actually  part  of  the  surface,  giving  rise  to  a 
phantom  surface  element.  In  principle,  this  should  not  occur  in  LSCM  and  one  might  expect  that 
it  would  be  possible  to  discriminate  between  real  surface  elements  and  phantoms  by  keeping  track 
of  the  intensity  associated  with  each  A(xi,  yj),  but  apparently  this  is  not  practical  and  frequently 
one  must  deal  with  images  containing  such  artifacts. 
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Since  the  artifacts  are  associated  with  geometrical  features  of  the  surface  being  imaged, 
they  are  approximately  reproducible.  Thus,  one  cannot,  in  general,  eliminate  them  from  an  image 
by  simple  statistical  techniques,  such  as  taking  a  median  elevation  from  a  number  of  images  of  the 
same  region  of  a  surface. 

Since  short-range  scaling  defines  the  fractal  dimension  of  a  fractal  surface,  one  cannot 
employ  smoothing  techniques  (e.g.,  low-pass  filtering)  to  deal  with  the  artifacts  on  fractal  surfaces 
without  significantly  altering  their  nature. 

Thus,  a  procedure  for  dealing  with  phantom  surface  artifacts  must  be  devised  in  order  to 
make  LSCM  a  practical  tool  for  quantitative  fractal  analysis. 

This  report  describes  two,  2-step  processes  for  the  treatment  of  phantom  surface  artifacts 
on  LSCM  images  of  fractal  surfaces.  Both  techniques  employ  "gradient  filtering"  and  "dilation" 
to  identify  parts  of  LSCM  images  that  are  free  of  artifacts.  One  technique  then  determines  the 
distribution  of  local  fractal  dimensions  only  in  parts  of  the  image  that  are  free  of  artifacts,  and  the 
other  "patches"  the  parts  of  the  image  that  contain  artifacts  before  performing  fractal  analysis. 

The  techniques  are  applied  to  LSCM  images  of  three  high-strength  and  high-toughness 
steel  alloy  Charpy  impact  fracture  surfaces  (ref  6),  one  fatigue  fracture  surface  of  the  same  steel 
alloy  and  one  ceramic  fracture  surface.  The  techniques  are  also  discussed  in  the  context  of 
mathematically  defined  Brown  surfaces  (refs  3,4)  and  detailed  results  are  presented  for  the  Charpy 
impact  fracture  surface  having  the  most  extensive  distribution  of  artifacts. 

EXPERIMENTAL  DETAILS 

The  Charpy  Test 

The  techniques  for  the  treatment  of  artifacts  on  LSCM  images  are  discussed  in  detail  for  a 
Charpy  impact  fracture  surface.  The  Chaipy  impact  test  (ref  6)  is  a  commonly  used  measure  of 
material  toughness.  In  a  Charpy  test,  a  square  v-notched  bar  is  struck  and  fractured  by  a  swinging 
arm  and  the  energy  absorbed  is  determined  from  the  amplitude  of  the  swing  of  the  arm. 

The  Example  Fracture  Sample 

The  example  Charpy  specimen  was  a  10-mm  square  v-notched  bar,  prepared  from  ASTM 
A723  (modified  AISI 4340  with  0.2  atomic  percent  vanadium)  steel,  which  had  been  heat  treated 
to  produce  tempered  martensite  having  nominal  strength  of  160  Ksi  and  hardness  of  38  on  the 
Rockwell  C  scale.  The  Charpy  toughness  of  the  example  fracture  surface  was  17.71.  The  other 
Charpy  fracture  surfaces  had  higher  toughness  values. 
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The  LSCM  Images 


All  of  the  LSCM  images  were  obtained  with  a  LaserTech  Confocal  Microscope  system. 
The  system  returned  a  476x640  array  of  integers  ranging  from  0  to  255  (gray  levels).  The  z- 
scales  were  selected  in  such  a  way  that  essentially  the  entire  range  of  256  gray  levels  was  required 
to  cover  the  imaged  surfaces  (including  artifacts). 

The  Example  LSCM  Image 

The  xy  spacings  were  1.0050  ^im  and  one  gray  level  difference  corresponded  to  1.2946 
|xm  for  the  example  fracture  sample.  Detailed  numerical  results  were  obtained  by  analysis  of  a 
257x257  region  in  the  LSCM  image,  which  was  selected  to  avoid  gross  surface  ^scontinuities. 
However,  essentially  the  same  (i.e.,  differences  in  <D>  were  less  than  0.02)  distributions  of  D- 
values  were  obtained  for  similar  procedures  on  arbitrary  257x257  regions. 

The  particular  40x40  pixel  region  of  interest  (ROI)  in  Figures  1  through  3  was  chosen 
because  of  the  prominent  feature  near  {i,j]  =  {110, 150}.  The  40x40  size  was  chosen  to  make 
the  features  reasonably  easy  to  see,  and  it  is  at  least  superficially  representative  of  40x40  sections 
on  the  640x476  LSCM  image. 

Figure  1  is  a  mesh  plot  (ref  7)  of  elevations  z(ij)  in  the  ROI  in  units  of  xy  spacings  (i.e., 
one  unit  corresponds  to  1.0050  |xm)  of  the  raw  LSCM  image  of  the  rough  (fractal)  fracture 
surface  contaminated  by  phantom  surface  artifacts. 

NUMERICAL  METHODS 

The  Triangulation  Fractal  Analysis  Algorithm 

A  detailed  description  of  the  triangulation  algorithm,  including  a  discussion  concerning  its 
sensitivity  to  uniformly  and  normally  distributed  random  noise,  is  presented  in  Reference  2.  The 
algorithm  determines  the  fractal  dimension  of  single-valued  surfaces  in  3-space.  The  idea  is  as 
follows: 

Consider  L  by  L  images,  where  L  is  taken  as 

L  =  2^  +  1  for  integer  n  (1) 

to  facilitate  triangular  approximations.  For  such  choices  of  L,  one  may  define  sets  of  triangles, 
which  tessellate  the  surface  with  no  gaps  or  overlaps,  whose  vertices  coincide  with  subsets  of 
points  on  the  lattice  and  whose  projections  in  the  xy-plane  are  right  triangles  of  side  lengths  Y{m) 
which  are  given  by 


Y{m)  =  2^  form  ^[0,1,.. .n] 


(2) 
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Then  the  small  Y(m)  approximations  to  the  fractal  surfaces  scale  according  to 


d  ln(A(m)) 
dln(r(m)) 


(3) 


where  Aim)  is  an  approximate  surface  area,  which  can  be  expressed  in  terms  of  individual  triangle 
areas 

for  the  triangulation  based  on  Yim)  by  Yim)  cells.  Equation  (3)  defines  the  fractal  dimension  D. 

A  slightly  more  general  form  of  the  triangulation  algorithm,  which  allows  for 
Y(a,b,c,...)  =  2a3b5^...fora  =  {0,1, ...nj},  b  =  {0,l,...n2},  c  =  (0,1  ....ns),...  (5) 


is  presented  in  Reference  2. 

The  Treatment  of  Phantom  Surface  Artifacts 

The  artifacts  are  treated  in  two  complementary  2-step  processes. 

Step  1.  Detection  of  Regions  Containing  Phantom  Surface  Artifacts 

The  regions  containing  artifacts  are  generally  delineated  by  the  following  sequence  of 

steps: 

1 .  Apply  a  two-dimensional  finite  impulse  response  filter  to  the  image.  In  particular, 
results  have  been  obtained  for  a  gradient  filter,  having  a  convolution  kernel  of  form 

a  1  1> 

1-8  1 

,111, 

and  for  Sobel  filters,  having  convolution  kernels  of  the  form 


f  1 

2 

n 

ri 

0  -1^ 

0 

0 

0 

or 

2 

0  -2 

^-1 

-2 

-1; 

a 

0  -1; 

Gradient  and  Sobel  filtering  are  standard  image  processing  tools  (ref  8). 
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2.  Choose  a  threshold  gradient  level  to  associate  with  artifacts  and  binarize  by  placing 
ones  where  the  absolute  value  of  the  gradient-filtered  image  exceeds  the  chosen  level  and  zeros 
elsewhere. 

3.  Expand  and  fill  in  the  high  gradient  regions  selected  previously  by  applying  one  or 
more  dilation  processes.  The  dilation  process  employed  adds  up  to  8-connected  pixels  (extra 
ones)  to  the  boundaries  of  the  objects  in  the  binarized  images.  We  refer  to  a  dilated  binarized 
image  as  a  mask.  Dilation  is  a  standard  image  processing  tool  (ref  8). 

Histograms  of  the  filtered-image  intensities  of  the  fracture  surfaces  studied  were 
continuous.  There  were  not  obvious  choices  for  the  threshold  levels.  Thus,  the  threshold  level 
and  the  number  of  dilations  were  selected  by  trial  and  error,  based  upon  visual  comparisons  of  the 
LSCM  image  and  the  mask.  A  satisfactory  mask  covers  all  those  pixels  that  contain  artifacts  and 
does  not  cover  an  excessive  number  of  uncontaminated  pixels. 

Satisfactory  masks  required  two  dilation  steps  and  covered  more  than  15  percent  of  the 
five  fracture  surface  images  studied  in  order  to  smoothly  cover  the  regions  that  contained 
artifacts.  The  artifacts  were  so  densely  distributed  on  the  example  Charpy  fracture  surface  that  in 
order  to  obtain  a  satisfactory  mask,  one  had  to  choose  a  threshold  level  that  yielded  a  binarized 
image  covering  10  percent  of  the  LSCM  image  before  dilation;  and  after  dilation  the  resulting 
mask  covered  34  percent  of  the  LSCM  image!  Exclusion  analysis  based  on  a  24  percent  mask 
yielded  <D>  larger  by  0.02  than  that  obtained  for  the  34  percent  mask.  Analysis  based  on  a  40 
percent  mask  yielded  the  same  <D>  to  three  digits. 

Figure  2a  is  a  MATLAB™  pseudocolor  plot  with  interpolation  shading  of  the  ROI  of 
Figure  1 .  The  gray  level  in  Figure  2a  is  proportional  to  the  measured  surface  elevation.  Figure  2b 
shows  a  contour  plot  delineating  the  34  percent  phantom  surface  artifacts  mask  for  the  same  ROI. 
The  34  percent  masks  obtained  by  the  gradient  and  the  two  Sobel  filters  were  in  essence  identical, 
and  employing  the  exclusion  or  the  patching  techniques  yielded  essentially  identical  statistical 
distributions  of  local  D-values.  The  correspondence  of  the  masked  regions  with  the  artifacts  may 
be  seen  by  comparing  Figure  2b  with  Figures  1  and  2a. 

Step  2a.  Exclusion  Analysis 

The  fractal  character  of  a  surface  could  be  described  exclusively  in  terms  of  local  D-values 
measured  outside  the  regions  containing  artifacts.  We  refer  to  such  a  procedure  as  exclusion 
analysis.  Exclusion  analysis  would  provide  a  complete  description  of  fi'actal  scaling  for  a  firactally 
homogeneous  surface.  We  refer  to  the  determination  of  the  distribution  of  local  D-values  based 
on  regions  entirely  outside  the  masked  portion  of  a  surface  under  study  as  "pure  exclusion 
analysis."  We  also  applied  a  form  of  exclusion  analysis,  which  we  refer  to  as  "modified  exclusion 
analysis,"  in  combination  with  the  patching  technique  described  below,  wherein  we  only  demand 
that  the  central  point  of  the  local  region  be  outside  of  the  masked  region. 
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Pure  exclusion  analysis  on  17x17  sections  excluded  62  percent  of  the  sections  for  the  34 
percent  mask  applied  to  the  example  Charpy  fracture  surface. 

Step  2b.  Local  Fill-In  Patching 

A  number  of  patching  schemes  were  evaluated  by  studying  masked  Brown  surfaces 
(refs  3,4).  The  local  fill-in  patching  technique  described  below  yielded  patched  surfaces  having  D 
within  ±0.04  of  the  D-value  measured  for  original  Brown  surfaces,  having  D  between  2.25  and 
2.65,  for  34  percent  masks  determined  for  the  example  Charpy  fracture  surface.  One  obtained 
closer  values  of  D  for  patched  Brown  surfaces  with  less  extensive  masks.  Lx)cal  fill-in  patching 
was  accomplished  by  an  iterative  procedure  in  which  the  masked  region  is  filled  in  by  following  a 
prescription  based  upon  locally  defined  elevations.  The  following  scheme  was  adopted: 

1.  Assign  "out-of-range  values"  to  those  pixels  in  the  raw  image  corresponding  to  ones  in 
the  mask. 

2.  Select  pixels  with  out-of-range  values. 

3.  If  at  least  5  of  the  8  nearest  neighbors  of  a  selected  pixel  have  "in-range"  values,  then 
assign  the  average  of  the  second  and  third  largest  near  neighbor  in-range  values  to  the  selected 
pixel.  If  less  than  5  of  the  8  nearest  neighbors  have  "in-range"  values,  go  to  the  next  selected 
pixel. 

4.  If  out-of-range  pixels  remain,  go  to  No.  2  above. 

Similar  procedures  were  employed  for  border  pixels.  For  example,  for  the  pixel  at  (i,/) 
define  the  8  nearest  neighbor  values  as  [z{l,i-l),  z{2,i),  z(2,i+l),  z(2,i-l), 

z{2,i+l),  z{2,i)},  and  patch  according  to  the  rules  for  interior  points.  However,  the  results 
reported  here  were  based  on  257x257  regions,  which  were  extracted  from  the  476x640  LSCM 
images,  and  were  selected  with  a  sufficient  number  of  bordering  pixels  to  eliminate  edge  effects. 

Essentially  the  same  distributions  of  local  £)-values  were  obtained  for  a  number  of  LSCM 
fracture  surface  images  and  Brown  surfaces  by  assigning  the  second  largest,  the  third  largest,  the 
average  of  the  second  to  fourth  largest,  or  the  average  of  the  third  and  fourth  largest  nearest 
neighbor  values. 

The  local  fill-in  patching  procedures  described  above  succeeded  in  defining  values  for  all 
masked  pixels  in  all  cases  studied.  Local  fill-in  patching  yielded  "reasonably  appearing"  patched 
fracture  surfaces  and  patched  Brown  surfaces.  Furthermore,  local  fill-in  patching  yielded  patched 
Brown  surfaces  exhibiting  approximately  the  same  fractal  scaling  as  the  original  Brown  surfaces. 
However,  the  34  percent  mask  patched  Brown  surfaces  were  visibly  different  from  the  original 
Brown  surfaces. 
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Figure  3  is  a  mesh  plot  of  the  ROI  after  the  2-step  treatment  based  on  the  34  percent  mask 
(ref  7).  Based  upon  the  results  for  patched  Brown  constructions,  the  patched  regions  of  Figure  3 
have  essentially  the  same  fractal  character  as  regions  outside  the  mask.  However,  the  patched 
regions  cannot  be  expected  to  be  faithful  representations  of  the  true  surface. 

APPLICATIONS  OF  EXCLUSION  AND  PATCHING  ANALYSIS 

Table  1  summarizes  the  results  for  the  example  Chaipy  fracture  surface  obtained  by 
analysis  of  the  raw  image  and  by  exclusion  and  patching  analyses.  The  values  of  <D>  obtained 
via  the  exclusion  and  patched  analyses  are  larger  than  the  coiresponding  "coastline  values"  Dc 
(viz.,  1.20  <  <  1.30)  obtained  by  slit-island  analysis  (ref  8)  of  Charpy  fracture  surfaces  of 

A723  steels  (ref  9).  However,  the  values  of  <D>  are  consistent  with  the  "coastline  values"  (viz., 
Dc  ~  1.40)  obtained  by  slit-island  analysis  of  fatigue  fracture  surfaces  of  A723  steels  (ref  9). 

Figure  4  exhibits  a  normalized  distribution  of  local  D-values  obtained  for  patched  surfaces. 
The  occupation  numbers  N  were  normalized  by  dividing  by  the  maximum  occupation  number. 

The  normalized  distribution  of  17x17  local  D-values  obtained  by  the  exclusion  procedure  (not 
shown)  is  essentially  the  same  as  that  obtained  for  the  17x17  patched  case. 

The  distributions  of  local  D-values  for  the  tougher  Charpy  fracture  surfaces  were  more 
complex;  The  17x17  exclusion-based  analyses  of  the  two  tougher  Charpy  samples  yielded  slightly 
larger  values  for  <D>  and  a  than  found  in  the  example  fracture  surface,  and  the  distributions  of 
local  D-values  exhibited  small  shoulders  on  the  low  D  side.  The  changes  in  <D>  and  a  were 
small  and  the  shoulders  were  subtle.  Analyses  of  the  patched  surfaces  were  consistent  with  the 
exclusion  results  in  the  17x17  cases.  Furthermore,  the  low  D  shoulders  became  more  distinct  in 
33x33  and  65x65  analyses  of  the  patched  surfaces.  However,  although  the  results  are  intriguing 
and  seem  to  suggest  a  change  in  mechanism  for  the  higher  toughness  Charpy  fracture  surfaces, 
one  must  remember  that  the  patched  regions  cannot  be  faithful  representations  of  the  actual 
fracture  surfaces. 

The  ceramic  and  the  fatigue  fracture  specimens  are  not  well  characterized.  However,  the 
following  results  were  obtained:  Exclusion  analysis  of  the  fatigue  and  ceramic  fracture  surfaces 
yielded  £>  =  2.37  ±  0. 10  and  D  =  2.28  ±  0.06,  respectively.  There  were  no  indications  of 
shoulders  on  the  distributions  of  local  D  values.  The  agreement  in  <D>  and  the  trends  in  o  were 
similar  to  those  in  Table  1. 

SUMMARY 

The  present  results,  in  combination  with  those  of  Reference  2,  support  the  conclusion  that 
quantitative  values  for  the  fractal  dimension  D  of  surfaces  may  be  obtained  by  analysis  of 
experimentally  acquired  images  contaminated  by  random  and  systematic  errors.  In  particular,  one 
can  obtain  useful  measures  of  D  from  LSCM  images  containing  substantial  phantom  surface 
artifacts. 
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The  exclusion  technique,  in  combination  with  the  triangulation  algorithm  or  a  similar 
fractal  analysis  algorithm  adapted  for  treatment  of  continuous  surfaces,  can  be  used  to  define  a 
distribution  of  local  D-values  characteristic  of  regions  of  LSCM  images  that  do  not  contain 
artifacts. 

The  patching  techniques  described  here,  in  combination  with  triangulation  or  another 
fractal  analysis  technique  that  is  adapted  to  the  treatment  of  continuous  surfaces  (e.g.,  slit-island 
technique  (ref  9)),  can  also  be  employed  to  obtain  useful  values  of  D.  The  patched  images  can 
also  serve  as  an  aid  in  visualization  of  surfaces,  particularly  when  the  patched  fraction  is  not  too 
large. 


Of  course,  if  the  surface  under  study  is  fractally  heterogeneous  and  the  regions 
contaminated  by  phantoms  are  different,  then  neither  exclusion  nor  patching  can  completely 
describe  its  fractal  properties. 
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Table  1.  Summary  of  Resulte  for  Example  Charpy  Fracture  Surface 


<D> 

a 

Local  Size 

S  urface/Procedure 

2.52 

257x257 

Raw  elevation  data 

2.52 

0.02 

129x129 

Raw  elevation  data 

2.52 

0.03 

65x65 

Raw  elevation  data 

2.52 

0.08 

17x17 

Raw  elevation  data 

2.39 

257x257 

Patched 

2.40 

0.01 

129x129 

Patched 

2.40 

0.03 

65x65 

Patched 

2.40 

0.07 

17x17 

Patched 

2.40 

0.07 

17x17 

Pure  exclusion 

2.40 

0.07 

17x17 

Modified  exclusion 

Note:  Values  of  <D>  differ  in  the  third  decimal  place. 
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Figure  2a.  A  pseudocolor  plot  of  raw  elevation  data  for  the  example  Charpy  fracture  surface 
ROI  of  Figures  1  through  3.  The  gray  level  is  proportional  to  the  elevation. 


Figure  2b.  The  34  percent  mask  for  the  ROI  of  Figures  1  through  3. 

The  masked  part  of  the  image  is  inside  the  closed  contours. 
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Normalized  N 


Figure  4.  The  distribution  of  local  D-values  for  the  patched  257x257  region  of  the 
"example  Charpy  fracture  surface."  The  distribution  of  local  D  for  17x17 
sections  is  essentially  the  same  as  that  obtained  by  exclusion  analysis. 
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